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Abstract: The real- and imaginary-time-formalisms of thermal field theory and 
their extension to gauge theories is reviewed. Questions of gauge (in-)dependence 
are discussed in detail, in particular the possible gauge dependences of the singu- 
larities of dressed propagators from which the quasiparticle spectrum is obtained. 
The existing results on next-to-leading order corrections to non-Abelian screening 
and dispersion laws of hard-thermal-loop quasiparticles are surveyed. Finally, the 
role of the asymptotic thermal masses in self-consistent approximations to thermo- 
dynamic potentials is described and it is shown how the problem of the apparently 
poor convergence of thermal perturbation theory might be overcome. 

1 Overview 

The theoretical framework for describing ultrarelativistically hot and dense 
matter is quantum field theory at finite temperature and density. At suf- 
ficiently high temperatures and densities, asymptotic freedom should make 
it possible to describe even the fundamental theory of strong interactions, 
quantum chromodynamics (QCD), through analytical and mostly perturba- 
tive means. These lectures try to cover both principal issues related to gauge 
freedom as well as specific problems of thermal perturbation theory in non- 
Abelian gauge theories. 

After a brief review of the imaginary- and real-time formalisms of thermal 
field theory, the latter is extended to gauge theories. Aspects of different treat- 
ments of Faddeev-Popov ghosts and different gauge choices are discussed for 
general non-Abelian gauge theories, both in the context of path integrals and 
in covariant operator quantization. The dependence of the formalism on the 
gauge-fixing parameters introduced in perturbation theory is investigated in 
detail. Only the partition function and expectation values of gauge-invariant 
observables are entirely gauge independent. Beyond those it is shown that 
the location of singularities of gauge and matter propagators, which define 
screening behaviour and dispersion laws of the corresponding quasi-particle 
excitations, are gauge independent when calculated systematically. 

At soft momentum scales it turns out to be necessary to reorganize per- 
turbation theory such that (at least) the contribution of the so-called hard 
thermal (dense) loops (HTL/HDL) is resummed. The latter form a gauge- 
invariant effective action, and their gauge-fixing independence is verified. The 
existing results of such resummations on the modification of the spectrum of 
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HTL quasi-particles at ncxt-to-leading order (NLO) are reviewed, and a few 
cases discussed in more detail, with special attention given to gauge depen- 
dence questions. 

It is shown how screening and damping phenomena become logarithmi- 
cally sensitive to the strictly nonperturbativc physics of the chromomagneto- 
static sector, with the exception of the case of zero 3- momentum. In particular 
the definition of a non-Abelian Debye mass is discussed at length, also with 
respect to recent lattice results. 

Real parts of the dispersion laws of quasiparticle excitations, on the other 
hand, are infrared-safe at NLO. However, as the additional collective modes 
of longitudinal plasmons and fermionic "plasminos" approach the light-cone, 
collincar singularities arise, signalling a qualitative change of the dispersion 
laws and requiring additional resummations. At high momenta only the reg- 
ular modes of the fermions and the spatially transverse ones of the gauge 
bosons remain and they tend to asymptotic mass hyperboloids. The NLO 
corrections to the asymptotic thermal masses play an interesting role in a 
self-consistent reformulation of thermodynamics in terms of weakly interact- 
ing quasi-particlcs. In this application, the problem of poor convergence of 
resummed thermal perturbation theory resurfaces, but it is shown that it 
may be overcome through approximately self-consistent gap equations. 



2 Basic Formulae 



Before coming to quantum field theories and gauge theories in particular, 
let us begin by recalling some relevant formulae from quantum statistical 
mechanics. 

We shall always consider the grand canonical ensemble, in which a sys- 
tem in equilibrium can exchange energy as well as particles with a reser- 
voir such that only mean values of energy and other conserved quantities 
(baryon/lepton number, charge, etc.) are prescribed through Lagrange mul- 
tipliers /3 = 1/T and on — —fi /Uj, respectively, where T is temperature and 
are the various chemical potentials associated with a set of commuting con- 
served observables iV; — Nj satisfying [Ni,Nj] = and [H,Ni\ = 0, where 
H is the Hamiltonian. 

The statistical density matrix is given by 



z- 1 



exp 



-aiNi - 0H 



z- 1 



exp 



-0(H - mNi) 
" v— — ' 

=: H 



(1) 



where Z is the partition function 

Z = Z(V,T, fH )=Tve^ B . (2) 

and the thermal expectation value (ensemble average) of an operator A is 
given by 

(A) = Tr[eA] . (3) 
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When total energy and particle numbers are extensive quantities^], i.e. 
proportional to the volume V, one also has In Z oc V, and since we shall be 
interested in the limit V — > oo, it is preferable to define intensive quantities. 
The most important one is the thermodynamic pressure 

p =w lnZ (4) 

which up to a sign is identical to the free energy density F/V = —P (F is 
also referred to as the thermodynamic potential S2). 

Other thermodynamic (or should one say thermo-static?) quantities can 
be derived from P, such as particle/charge densities 

1 - dP 

rn = v (N i} = -, (5) 
energy density 

-*/"4<*>-4^~T 

(at fixed ai), and entropy density (which will play a prominent role at the 
very end of these lectures) 

a = S/V = ~<- Ing) = - In Z + ^{H - ^N t ) 
dP 

= ^=P{P + e- fum) (7) 

In the latter equations one recognizes the familiar Gibbs-Duhem relation 

E = -PV + TS + fiiNi, (8) 

which explains why P was defined as the (thermodynamic) pressure. A pri- 
ori, the hydrodynamic pressure, which is defined through the spatial compo- 
nents of the energy-momentum tensor through |(T"), is a separate object. In 
equilibrium, it can be identified with the thermodynamic one through scal- 
ing arguments |9^| , which however do not allow for the possibility of scale 
(or "trace" ) anomalies that occur in all quantum field theories with non-zero 
j3- function (such as QCD). In |5l| it has been shown recently that the very 
presence of the trace anomaly can be used to prove the equivalence of the 
two in equilibrium. 

All the above relations continue to hold in (special) relativistic situations, 
namely within the particular inertial frame in which the heat bath is at 
rest. In other inertial frames one has the additional quantity of heat-bath 
4- velocity and one can generalize the above formulae by replacing V = 
J d 3 x —>■ J dE^u^, and the operators in (|l]) by 



dS^Uv, AW dEfd?. (9) 



1 A notable exception occurs when general relativity has to be included. 
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a.i and [3 are Lorentz scalars (i.e., temperature is by definition measured in 
the rest frame of the heat bath). The partition function can then be written 
in covariant fashion as |69l 



Z = Tr 



exp IdEA-T^Pu-ti 



(10) 



where we have introduced an inverse-temperature 4- vector = 

However, in what follows we shall most of the time remain in the rest 
frame of the heat bath where = 5^. 



3 Complex Time Paths 

With a complete set of states given by the eigenstates of an operator tp, 
<p\<p) = ip\tp), we formally write 



Z = Tr 



e -/MT 



= ^(<p\e-f> A \<p). (11) 



In field theory, we shall boldly use the field operator tp = tp(t, x) in the 
Heisenberg picture and write \tp) for its eigenstates at a particular time. 

Using that the transition amplitude (tpi |e _lff ( il_t °) \tp ) [the overlap of 
states that have eigenvalue <po(x) at time t with states that have eigenvalue 
tpx(x) at time ti] has the path integral representation 



(tpt\e- iH ^-^\tp 



to 



M / Dtp expi / dt / <TxC{tp,dtp) (12) 



we can give a path integral representation for Z that takes care of the density 
operator by setting t\ — to — i/3, and of the trace by integrating over all 
configurations with tpx = tp . 

When there is a chemical potential /i ^ 0, we have H — > H = H — (iN 
and this implies L — > L + fiN if there are no time derivatives in N, and we 
have 

//■to— ifi r 
Vtp expi J dt J d 3 x£ (13) 

where the path integral is over all configurations periodic in imaginary time, 
<p(t ,x) = tp{t - if3,x). 

In this formula, real time has apparently been fixed to to and replaced 
by an imaginary time flow which is periodic with period /3, the inverse tem- 
perature. In equilibrium thermodynamics, this seems only fitting as nothing 
depends on time in strict equilibrium. 
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However, we have not really required time to have a fixed real part. We 
have made the end point complex with the same real part, but the integration 
over t in need not be a straight line with fixed to . Instead we shall consider 
a general complex time path, and this allows us to define Green functions by 
the path integral formula 

(T c ifi ■ ■ ■ <p n ) — A/" J T>ip <pi ■ ■ ■ (p n exp i J dt J d 3 x C (14) 

c 

where T c now means contour ordering along the complex time path C from 
to to to — i/3 such that ti £ C, and t\ y t% y ■ ■ ■ y t n with respect to a 
monotonically increasing contour parameter. 

Through quantities like ( |l4| ) we are not restricted to time- independent, 
thermo-static questions, but may also consider small perturbations of the 
equilibrium (response theory). 

The time path introduced in (14) is in fact not completely arbitrary. 



Considering spectral representations in the energy representation leads to 
the conclusion (9^] that C has to be such that the imaginary part of t is 
monotonically decreasing. This is a necessary condition for analyticity; in 
the limiting case of a constant imaginary part along (parts of) the contour, 
distributional quantities (generalized functions) arise. 

Except for the periodic boundary conditions with regard to the end points 
of C (which become anti-periodic for fermionic field operators and Grassmann- 
valued "classical" fields), the path integral formula (|lj) is formally identical 
to the familiar one from T — and /i = 0. 

Indeed, perturbation theory is set up in the usual fashion. Using the 
interaction-picture representation one can derive 

(TA--0„)-|(TA--One i ^ I )o, (15) 

where Ci is the interaction part of £, and the correlators on the right-hand- 
side can be evaluated by a Wick(-Bloch-DeDominicis) theorem: 

(T c e i /c d4 ^^) =exp{-i / f d A x d 4 x'j{x)D c (x - x')j(x')}, (16) 
2 Jc Jc 

where D c is the 2-point function and this is the only building block of Feyn- 
man graphs with an explicit T and [i dependence. It satisfies the KMS (Kubo- 
Martin-Schwinger) condition 

D c {t - i/3) = ±e-^D c (t), (17) 

stating that e lflt D(t) is periodic (antiperiodic) for bosons (fermions). 
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3.1 Imaginary-Time (Matsubara) Formalism 

The simplest possibility for choosing the complex time path is the straight line 
from to to to — i/3, which is named after Matsubara ]9j who first formulated 
perturbation theory based on this contour. It is also referred to as imaginary- 
time formalism (ITF), because for to — one is exclusively dealing with 
imaginary times. 

Because of the (quasi-)periodicity (|l7|), the propagator is given by a 
Fourier series 

D%t) = ^-y j D(z u )e- iz " t , D(z v )= [ 1 dtD c (t)e iz ^ (18) 
-IP " Jo 

with discrete complex (Matsubara) frequencies 

_ . , n fZ bos. 

^ = W/3 + M, ^ e (z-iferm. (19) 

The transition to Fourier space turns the integrands of Feynman diagrams 
from convolutions to products as usually, with the difference that there is no 
longer an integral but a discrete sum over the frequencies, and compared to 
standard momentum-space Feynman rules one has 

/ - r 1 E / -§^, i(2^ 4 w - m^sHk). ( 20) 

However, all Green functions that one can calculate in this formalism are 
initially defined only for times on C, so all time arguments have the same real 
part. The analytic continuation to different times on the real axis is, however, 
frequently a highly non-trivial task p2| , so that it can be advantageous to 
use a formalism that supports real time arguments from the start. 



3.2 Real-Time (Schwinger-Keldysh) Formalism 

In the so-called real-time formalisms, the complex time path C is chosen such 
as to include the real-time axis from an initial time to to a final time tf. Since 
we have to end up at to — i/3, this requires a further part of the contour to 



run backward in real time |11S,19| and to finally pick up the imaginary time 
—i/3. There are a couple of paths C that have been proposed in the literature. 
The oldest one due to Keldysh J78| is shown in Fig. and this is also (again) 
the most popular one. 

Clearly, if none of the field operators in (|l4|) has time argument on C± or 
C2, the contributions from these parts of the contour simply cancel and one 
is back to the ITF. 

On the other hand, if to — ► —00, and all operators have finite real time 
arguments, the contribution from contour C3 decouples because from the 



Thermal Gauge Field Theories 7 

Im t 



<— to 


Ci 

>■ 


tf —* +00 










c 3 


to - i/3 


c 2 


Re t 



Fig. 1. Complex time path in the Schwinger-Keldysh real-time formalism 



spectral representation one has for instance for the propagator connecting 
contour C\ and C3 

00 

D^(k,t-(t -i\)) = J ^e- i ^)_^_ 0(fcjW) ^oo o (21) 

— OO 

for A G (0,(3) by Riemann-Lebesgue 

With only C\ and Ci contributing, the action in the path integral decom- 
poses according to 

/>oo />oo 

C(<p)= dtC{v (1) )- dtC{^ 2) ) (22) 



C1UC2 

where we have to distinguish between fields of type 1 (those from contour C\) 
and of type 2 (those from contour C2) because of the prescription of contour 
ordering in (|l4|)P| From (|2^) it follows that type-1 fields have vertices only 
among themselves, and the same holds true for the type- 2 fields. However, 
the two types of fields are coupled through the propagator, which is a 2 x 2 
matrix with non- vanishing off-diagonal elements: 

, _ ( (T V (i) V (t<)) a(cp(t')<p(t)) \ 

Here T denotes anti-timeordering for the 2-2 propagator and a is a sign which 
is positive for bosons and negative for fermions. The off-diagonal elements do 
not need a time-ordering symbol because type-2 is by definition always later 
(on the contour) than type-1. 

2 There are cases where this line of reasoning breaks down. The decoupling of 
the vertical part of the contour in RTF does however take place provided the 
statistical distribution function in the free RTF propagator defined below in 

§does have as its argument |fco| and not the seemingly equivalent \fk 2 + m 2 

^ Type-2 fields are sometimes called "thermal ghosts" , which misleadingly suggests 
that type-1 fields are physical and type-2 fields unphysical. In fact, they differ 
only with respect to the time-ordering prescriptions they give rise to. 
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In particular, for a massive scalar field one obtains 

where <5 ± (fc 2 — m 2 ) = 9(±ko)S(k 2 — m 2 ) . The specifically thermal contribution 
is that of the second line. Mathematically, it is a homogeneous Green function, 
as it should be, because it is proportional to S(k 2 — m 2 ). Physically, this part 
corresponds to Bose-Einstein-distributed, real particles on mass-shell. 

The matrix propagator ( ^4|) can also be written in a diagonalized form 

is, 

D c (fc)=M(fc )(^ JL)M(fc ) (25) 
with G F = l/(k 2 + m 2 + ie) and 

In the T — > limit one has 

M W ^Mo(fco)=(^ o) (27) 

so that one still has propagators connecting fields of different type. However, 
if all the external lines of a diagram are of the same type, then also all the 
internal lines are, because Yii = when k9^ = and any connected 

region of the other field-type leads to a factor of zero. 



4 Gauge Theories — Feynman Rules 

As a simple application of the formalism developed above and as a demon- 
stration of the need for more formalism for gauge theories, let us try to 
calculate the thermodynamic pressure of photons in the imaginary-time for- 
malism and in a covariant gauge. (There is no need for the real-time formal- 
ism here, because we are not considering Green functions external lines.) 
The simplest gauge to perform calculations is usually Feynman's gauge, 
which simplifies the Lagrangian of the electromagnetic fields according to 
C = -\F^F^ — \A V UA". 

One would therefore expect the partition function to be given by 

-i/3 

I VA cxpi J dtd 3 x£ = const. x(dctD)pj i ^ ic , 

J periodic J 
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and the thermal pressure would be calculated from 
In Z = -4 x ^Tr In □ + const 



= -4 x \vY, J ^3 ln K + fc2 ) + const 

J (27T 



Pk ln(l-e"^) 



const 



2 

as 

P(D - P(0) = ±lnZ = 4Tj^- 3 H l - e-^ = 2 x ^ 

giving twice the correct result of Planck for blackbody radiation. 

The error we made is that we have not calculated Tre - ^ in a physical 
Hilbcrt space (in fact in no Hilbert space at all, because there are negative- 
norm states). Instead of only two physical (transverse) degrees of freedom, 
we have added up the contributions from four. The standard way to get rid 
of the unphysical degrees of freedom is to cancel their contribution by ghost 
contributions, which evidently are required already in the Abelian case. 



4.1 Path Integral Faddeev- Popov Trick 

Because of gauge invariance under 6A^ = D" l b (A)ui b (where a, b are possible 
color indices), there is a redundancy in the path integral that leads to zero 
modes in the kinetic kernels, making them non-invertible. Thus, in order to 
be able to write down propagators and do perturbation theory, one needs 
to remove this redundancy. Using a suitable gauge fixing function F a [A](:r), 
this can be done by inserting 

r) F a 

n^MO^-det^ (28) 

into the measure of the path integral, selecting only one gauge field configu- 
ration per gauge orbit. 

One can equally well use F a [A](x)—( a (x) in place of F a [A](x) with arbi- 
trary functions ( a and perform a Gaussian average 

J VQe^ 2 ... 

over the latter. This gives so-called general or inhomogeneous gauge breaking 
terms 

£ —> C + —F a [A] 2 (29) 
2a 

with gauge fixing parameter a. 
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In covariant gauges F a [A](x) — d^A^(x) and Abelian electromagnetism, 
where a takes only one value and D^ b (A) — > d^, the determinant in (p8| ) is 

9F i 
det — = det □ = (det □)+ 2 x2 

so this indeed compensates for the two unphysical degrees of freedom in the 
above miscalculation of blackbody radiation. 

Usually, this "Faddeev-Popov determinant" does not play a role in QED 
because it is field independent. For calculating thermodynamic potentials, it 
does, because this determinant depends on temperature through the bound- 
ary conditions. 

In non- Abelian gauge theories, det |^ = det (^ F D b ^(A)j is field de- 
pendent, and it is convenient to introduce anticommuting and realQ Faddeev- 
Popov ghost fields 

VcVc expi J c c a ^-D b ;(A)c c = const, x det ^—D^A)) . (30) 

The correct boundary conditions are clearly those of the gauge potentials 
and thus are periodic in imaginary time despite the fact that the ghosts are 
anti-commuting and thus behave like fermions with regard to combinatorial 
factors in front of Feynman diagrams p(| . 



4.2 Covariant Operator Quantization 

While the path integral makes it evident how to treat ghosts at finite tem- 
perature, one can arrive at the same conclusion without recourse to path 
integrals in covariant (BRS) operator quantization |36|]87|, where at first it is 
somewhat surprising that anticommuting fields should end up having periodic 
rather than antiperiodic boundary conditions in imaginary time. 

BRS quantization is preferably done with Lagrange multiplier fields B 
and a gauge-fixed Lagrangian (in general covariant gauges) 

C = C inv - A a ^B a + ^B a B a - i(<^c a )ZV (31) 

with c and c being anticommuting field operators. 

The gauge-fixed Lagrangian possesses a global fermionic (BRS) symmetry, 
which in any gauge reads 

[iQbrS; Afj] = D^c, {IQbrs, c} = -§ c x c, , . 

[iQ B RS,£] = 0, {iQbrs, c} = iB, 

where we used a vectorial notation to write for instance D^c = (d^+gA^ x )c. 
In covariant gauges for instance this is generated by 



?brs = d x 



B ■ D a c - c ■ d B + ^(d c) • (c x c) 



(33) 



c is not the conjugate of c, but an independent field. 
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The BRS operator is nilpotent, Q BRS = and commutes with Lagrangian 
and Hamiltonian, [iQbrs,£} = = [iQbrs,^]- 

There is one further global symmetry, ghost number, with conserved 
charge 

N c = J <fx [doc-c-c- D a c] (34) 
satisfying 

[N c ,c]=c, [N c ,c]=-c, [JV c ,^]=0, [N c ,B]=0. (35) 

N c is anti-hermitian, N c = —NJ, although it has real eigenvalues n g h, 
which is possible because our arena is a non-Hilbert space V containing 
negative-norm states. 

The negative-norm states can be eliminated by the linear condition 

V -> Vphy S : QBRslphys) = 0, (36) 

and the true physical Hilbert space is finally obtained by modding out zero- 
norm states, 



~Hphys = Vphys/Vo . (37) 

The corresponding projection operator V in H p \ lys = VV can be shown 
[ p6p7[ | to have a complement that is "BRS exact" , meaning 

V + {Q BRS ,TZ} = 1, (38) 

but the actual construction of these operators is rather complicated. 

However, we apparently need them to be able to define the trace restricted 
to the physical Hilbert space that appears in Z = Trl e~@ H = TrT'e - ' 3 ^ 
or in expectation values of observables. 



Hata-Kugo Trick There is however an elegant trick that avoids the explicit 
construction of V Q: From [N c , Qbrs] = Qbrs it follows that A^Qbrs = 
Qbrs(^c + 1)" and therefore e^ N ^Q BRS = Q B Rse i7r(JV ° +1) = -Qbrs^, so 

{e^,Q B Rs}=0. (39) 

This, together with N c \^p) — for \ip) 6 H p hy S can be used to write 

Z = TrVe- fm = TrV e mNa e~ pH 

= Tre inN °e- 0H - Tr{Q B RS,^} e^e""* (40) 

Tr^{e i7rJV %Q B R S }e-^ 

S v ' 

o 

since [<2brs,#] = 0. 
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So the comparatively simple operator e 17rJVc can be used in place of the 
complicated V to express Z, and similarly thermal expectation values of 
gauge-invariant observables 0, through a trace in unrestricted V containing 
ghosts and other unphysical degrees of freedom, 

Z = Tr [^ N 'e- fiS ] , (O) = Z^Tt [e mN ^- 0H O] . (41) 

This result shows that the anticommuting ghosts, which naturally are sub- 
ject to antiperiodic boundary conditions, acquire a purely imaginary chemical 
potential 

tic = in/p. (42) 
In the ITF, the Matsubara frequencies of the ghosts are therefore 

z v = 2vri(ri - -)/j3 + [i c = 2irin/f3, neZ (43) 

like those of ordinary bosons. Thus, while they do have fermionic combina- 
torics in Wick contractions and the like, thermodynamically they behave like 
bosons. 

In the RTF, where the matrix-valued propagator ( p4| ) can be written as 

-iD=( G * _° G *)+(G F -G* F )x 

f 6(k )n(kq) + 9{-k )n(-k ) s_gn(fc )n(fc ) \ (aa) 

X<7 { sgn(fc )(a + n(fc )) 0(k o )n(k o ) + 9(-k Q )n(-k ) J [ * V 

with a — ± for bosons/fermions, n(x) — , and fco — — Mi we have 

a = — for the Faddeev-Popov ghosts, but 

n FT >(ko - m/p) = ePko } i7T - 1 = -n B E(k ), (45) 



so the imaginary chemical potential ( |42[ ) in effect negates a and replaces 
Fermi-Dirac by Bose-Einstein statistics. 

Only the signs arising in Wick contractions are those of fermions, which 
shows that the Faddeev-Popov ghost propagators have indeed the right form 
to be able to compensate for unphysical degrees of freedom contained in 
the gauge boson propagator, which naturally has a = + and Bose-Einstein 
statistical factors. 

Compared to (H), the gauge boson propagator also has a factor 9^ = 
{~9iiv + (1 — a) -rs^ ) in covariant gauges. We shall also consider other gauges 
in what follows, which can be characterized by 

g^^5^=g^ ~— + (/ -ak ) (46) 

where / is the momentum-space version of / in .F a [A] = f^A" 
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Popular gauge choices besides the familiar covariant gauges include axial 
gauges (F a — n^A®, const.) and Coulomb gauge(s) (F a = d l A"). 

In axial gauges, ghosts decouple completely, because the Faddeev-Popov 
determinant det(n-d) is field- and temperature-independent, however they are 
fraught with technical difficulties, already at zero temperature. The particu- 
larly attractive "temporal" gauge m) 1 = Sq, which does not cause additional 
breaking of Lorentz symmetry, is unfortunately inconsistent with periodic 
boundary conditions. Relaxing those leads to rather complicated Feynman 
rules in the ITF jn|, while the RTF version seems more tractable |fro| , at 
least it does not appear to be more problematic than at zero temperature. 

Coulomb gauge is in fact widely used at finite temperature, because it also 
does not lead to additional Lorentz symmetry breaking. However, it does have 
less simple Slavnov- Taylor identities Q because ghosts do not decouple, 
although they frequently do not contribute, since their (RTF) propagator 
does not contain statistical distribution functions. 

4.3 Frozen Ghosts 

In J90| , alternative Feynman rules have been proposed which avoid thermal- 
ized ghosts even in covariant gauges. In the RTF one can switch off the 
interactions as to — > — oo, and define the physical Hilbert space in terms of 
Abelianized in-states. Without additional Lorentz symmetry breaking, phys- 
ical in-states can then be chosen as 

|phys, in) = \A phys . -quanta) |0 w.r.t. A unphys .,B, c, c) (47) 
with ^phys.( fc ) =W(k)A«(k) and 

A " = 0, A 11 = -(S ij - (48) 

The unphysical states correspond to 

^,.(i) = (f-r(iK(*), B\ c°,andc». 

The interaction-picture free Hamiltonian then separates into two com- 
muting parts 

Hqi = # P / ys + # 7 phyS , [#o/ ys , Hfi phy *\ = 0, (49) 
and, because 7JQ" phys |phys, in) = 0, thermal averages factorize: 

2J (phys, in|e _/3HoJ • • • A phys . ■ ■ ■ A unphys . • • • c • • • c • • • |phys, in) = 

phys 

^2(A phys .\e- f3H ^ yi " ■ ■ ■ A phys . ■ ■ ■ |A phys .}(0|- --Amphys. ■■■C---C-- -|0) 

phys 
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with the thermal Wick theorem applying to the first factor under the latter 
sum, and the T — Wick theorem to the second one. 

This leads to alternative Feynman rules for gauge theories in RTF in which 
only the transverse projection of the gauge bosons have thermal (matrix) 
propagators, and all other fields have only the T = limits of those. The rest 
of the Feynman rules is as usual in RTF. 

E.g. in Feynman gauge (a = 1) the gauge boson and ghost propagators 
now read 

D^ = -i^M^ G Q F _y M-iG^-^Mo ( G Q F _g*JMo 
= -.^( 2 ™ 2) 2nS y 2) y2n6(ki)A»M\ko\)({ J), (50) 

In general linear gauges one has to replace in the vacuum part according 
to (fi"6|); the thermal part remains unchanged. 

Using these Feynman rules simplifies certain calculations in covariant and 
other gauges Q, because, although ghosts are present, they do not carry sta- 
tistical distribution functions, but are "frozen" . On the other hand, the usual 
cancellation of pinch singularities in the RTF (absence of "pathologies"), 
turns out be more complicated, and occurs in general only upon Dyson re- 
summations Pll. 



5 Gauge Dependence Identities 

As we have seen, perturbation theory and its Feynman rules require the 
introduction of gauge fixing terms. Clearly, physical results have to come 
out independent of those. We shall therefore now study in detail to what 
extent perturbative calculations will exhibit dependences on the gauge fixing 
parameters by considering 

F a [A] — > F a [A] + SF a [A] . (52) 

With <5F°L4] tx F a L4], this also comprises a possible variation of the gauge 
parameter a in (|9|) or (|3l"|). 



5.1 Gauge Independence of the Partition Function 

Path Integral The introduction of the gauge fixing term together with 
the Faddeev-Popov determinant according to ( p8| ) was done in a way that 
picks one representative field configuration from each gauge orbit.^j So by 

5 At least perturbatively this is guaranteed by the existence of dF/dui; non- 
perturbatively there may be obstructions to worry about. 
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construction the partition function or averages of gauge-invariant operators 
are independent of the gauge fixing terms. 

This can be checked explicitly by noting that the variation (|5^) of the 
gauge fixing function can be written as 

-l 

SF[A}. (53) 



9F W D[A] 



dA 



k[A] 

A corresponding change of the gauge breaking term -^{F a ) 2 is thus equiv- 
alent to a gauge transformation SA^ — Z3 M [v4]<5£ with the above non-local, 
field-dependent parameter <5£ [A]. 

The invariant part of the action is of course invariant under A — > A+5A, as 
are any gauge invariant operators that might have been inserted, so it remains 
to check that the path integral measure T>A together with the Faddeev-Popov 
determinant is invariant, too. This can indeed be verified by writing ST>A as 
trdSA/dA, and 5det[ff • D] = det[|f • D] x £(trlog[§§ • D]) = det[ff • D] x 

* r ([f5 ' ^lZ. "I§3 ' an< ^ then using that gauge transformations form a 
group MM. 



Covariant Operator Formalism In the covariant operator formalism the 
Lagrangian in a general gauge F can be written as 

£ = £ hw + B-F+^B-B-c-[Q BRS ,F] (54) 

and variations of F correspond to BRS transformations with parameter c-5F: 

6£ = B-5F-c- [Qbrs , SF] = {Qbrs ,c-6F} (55) 

It is always possible to construct an operator 8G such that SH = {Qbrs, SG} 
(if no time derivatives are involved, one simply has 6G = — Jd 3 xc ■ SF). 

Gauge independence of the partition function and of thermal averages of 
gauge invariant operators defined by ( |4l| ) can be verified as follows 66 : 

Variations of exponentiated operators can be expressed as 

A 

e A+SB _ e A = J dx e XA SB e (l-X)A + 0( £ 2) (56) 


and using this one finds 

STr[e- pB+iKN '0] = - / d\ Tr[c~ Aff {Qbrs, SG}e XH e~^ H e i7TNc O] 
Jo 

= -/ d\Tr[e- XH 5Ge XH {e iwA S Q BRS } e"^©] (57) 

Jo s v ' 



because of cyclic invariance of the trace and [Qbrs >0] — for a gauge- 
invariant operator O. 
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5.2 Gauge Dependence of Green Functions 

Green functions, i.e. thermal averages of products of field operators like 
(T c • • • ■ ■ ■ tp ■ ■ ■ if) ■ ■ •), are, however, gauge-variant objects, and will there- 
fore in general depend on gauge fixing parameters.^ 

In particular, the propagators of gauge and matter fields will contain all 
sorts of gauge parameter dependences. Yet, they are among the prime objects 
of linear response theory as they are used to derive the properties of quasi- 
particles. 

Historically, a stimulating failure was the attempt to extract the damp- 
ing constant of long-wavelength plasmons in a gluon plasma from one-loop 
thermal perturbation theory. Some of the results that were accumulated in 
the 80 's are summarized in Table |l|. These turned out to be gauge indepen- 
dent in algebraic gauges but gauge-parameter dependent in covariant ones. 
Moreover, in the latter the damping constant came out with the wrong sign 
which some took as signal of an instability of the perturbative ground state 




Table 1. Bare one-loop gluonic plasmon damping constant 7(|fc| — ► 0) 



gauge 



published 



-[T + (§ -2) 2 ] 
+ 1 
+ 1 

[11 + 1(1 -a) 2 ] 

_45 
4 

-11 



covariant gauges 
temporal gauge 
Coulomb gauge 
background covariant gauges [ ^ 



gauge-independent effective action [ 1.22,112 
gauge-independent pinch technique Q 



1980 [75] 
1985 72 1 
1987 [67] 
1987 [65] 
1987 [65] 
1988 [F02] 



It was in particular Pisarski [108] who argued that these results were 



simply incomplete, and who together with Braaten |3^,|39|,|36| devised an ap- 
propriate resummation scheme. However, since explicit calculations can only 
be performed in practice with a rather limited choice of gauge parameters, 
it is important to investigate more generally how gauge dependent the full 
propagators are and whether they contain gauge-independent information at 
all. To this end, we shall first derive rather general "gauge dependence identi- 

6 Notice that gauge invariance and gauge-fixing parameter independence are sep- 
arate issues: a functional of fields can be gauge invariant and yet depend para- 
metrically on the gauge-fixing function; conversely, independence of gauge fixing 
parameters (within a class of gauges) does not imply that a particular functional 
(e.g. of mean fields) is a gauge-invariant one. 
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a 1 

j<1 —a D a [<Pl 5~0 8F a [ip] 

Fig. 2. Additional Feynman rules for <5X J [<y3] in the case of linear gauge generators 
D and linear gauge fixing F 



ties" and study their consequences for the thermal Green functions of interest 



Primary Diagrams In order to unclutter the relevant relations, we shall 
temporarily switch to the compact notation of DeWitt ^7j, where a single 
index i comprises all discrete and continuous field labels (e.g. i = {A, /i, a, x)) 
and Einstein's summation convention is extended to include integration over 
all of space and time (the latter along the contour C). This way, ip 1 will 
represent an arbitrary gauge or matter field ip 1 — {A®* (x^, ip®* (xi), . . .}; only 
the Faddeev-Popov ghosts fields will be treated separately. 
The generating functional of Green functions reads 

Z[J] = (e iJ ^) with Jap* = Jd 4 x[J {A}tia (x)A^ a (x) + ...} (58) 

c 

and depends implicitly on a gauge fixing functional F a [ip], where a — (a Q , x a ) 
comprises both a group and a space-time index. 

Information on the dependence on F a can be obtained either by using 
BRS techniques or equivalently by employing the non-local gauge transfor- 
mation of (^3|), which in compact notation reads 

Sip* = D^MSF with 6C = 5t a [<p] = -9 a f l [<p]6Ff , [(p] (59) 

where D l a is a generalized function containing the gauge generators, and 
S a /3M = ~{FiD i ol )~ 1 [<p] is the Faddeev-Popov ghost propagator in a back- 
ground field ip. This immediately gives 

8\nZ[J] =U j (Dl t [ip]T[3Vp\5F () [ip])[J] under F a — > F a + 5F a . (60) 

The diagrammatic content of ( |60| ) is more conveniently investigated after 
a Legendre transformation of W[J] = —ilnZ[J] to the effective action 



6W[J] 



r[<p]=W[J]-J j <pi, = (61) 



which is the generating functional of one-particle-irreducible (1-p-i) diagrams. 
Equation (|60|) then becomes 



Sr[0[ = ^1 {DiMTM SF p [tp]) m = rM 6Xi[Cp] . (62) 
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1^ 

+ .<-^----^) + i - <-^— --^) + • • ■ 

Fig. 3. Primary diagram expansion of 5X 3 [tp] through 2-loop order. Contributions 
involving an undifferentiated 8F a [<p] have been dropped, which corresponds to omit- 
ting the trivial tree-level gauge dependence [ p3[ 

Diagrammatically, Jj [0\ is the sum of all mean-field dependent (primary) 
1-p-i one-point diagrams, while JA^i/j] is given by primary diagrams which 
involve the additional vertices introduced in Fig. ||and which are 1-p-i except 
for the basic ghost line attached to SF 13 , as shown in Fig. 

From these relations one can derive gauge dependence identities for 1-p-i 
vertex functions by differentiation with respect to (p. For example, the gauge 
dependences of the 2-point vertex function (self-energy) r^j [0] are determined 
by the diagrams shown in Fig. ^. 



QED As a first application let us consider the case of an Abelian theory such 
as QED. The additional Feynman rules of Fig. ^ involve the ghost propagator, 
but there are no further ghost vertices in the theory (for linear gauge fixing) , 
so only the very first diagram in Fig. || arises. 

Furthermore, the structure of the gauge generator is such that 

'«—{?• lr A " w ./-={ ±ie tf4iJ „v <«> 

If the external indices i,j of ST^j correspond to photons, one finds that 
one cannot even build the one remaining diagram of Fig. 0, so STa^a* proves 
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Fig. 4. Gauge dependence ol the 2-point vertex function Aj'[0] assuming that all 
one-point functions vanish at (p — 

to be completely gauge-fixing independent. This is in fact a well-known result 
which can be understood also by the gauge invariance of the electromagnetic 
current operator. 

On the other hand, if the external lines are fermionic, there is a non-trivial 
right-hand-side to Fig. |], as shown in Fig. [|, so the fermion self-energy is a 
gauge fixing dependent quantity, already in QED. 




Fig. 5. Gauge dependence identity for the fermion self-energy 



Hard Thermal/Dense Loops In the high-temperature (large-chemical- 
potential) limit of QED and QCD, it turns out that the leading contributions 
to the 1-loop vertex functions, the so-called "hard thermal (dense) loops" 
(HTL/HDL) obey tree-level- type ghost-free Ward identities and appear to 
be gauge-fixing independent ^9|j3|^]. This gauge independence, however, 
does not arise in an obvious way and involves non-trivial cancellations in the 
various gauges that have been considered. 
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The above gauge dependence identities can be used to verify the gauge 
independence of the HTL's in a rather simple manner. The only further 
ingredients needed are the temperature power-counting rules given in J39j, 
which, roughly, read as follows: in a Feynman diagram, explicit loop momenta 
in the numerator give a factor T, each propagator counts as T -1 , and the 
sum-integral over the loop momentum contributes T 3 unless there are two 
or more propagators with the same statistics, in which case the sum-integral 
counts as T 2 . 

By this, the leading temperature contributions to a 1-loop vertex func- 
tion are found to be proportional to T 2 , such that an iV-point gluon vertex 
function scales as fnv) ~ g N T 2 k 2 ~ N (where k represents generically com- 
ponents of external momenta). If two external lines are fermionic, we have 
r,(N) ~ g N T 2 k 1-N , while vertex functions with more than two external 
fermion lines do not contribute terms oc T 2 . 



Considering e.g. vertex functions with only external gluons, all of the 
potential gauge dependences of the HTL's are contained in the 1-loop contri- 

N 




The above temperature power-counting rules are modified only by the ad- 
ditional vertex 8F" [0] which may bring in one power of loop momentum and 
thus one power of T in derivative gauges, or none in algebraic ones. Adding 
up, one finds that the right-hand-side of (p4|) is proportional to T° - 1 . Hence, 
gauge dependences of one-loop vertex functions can occur only at sublead- 
ing order cx T. HTL(HDL)'s on the other hand are found to be completely 
gauge- fixing independent. 

5.3 Gauge Independence of Propagator Singularities 

In non-Abelian gauge theories, all the matter and gauge field vertex func- 
tions, self-energies, and propagators contain highly nontrivial gauge depen- 
dences, which raises the question whether there is any gauge-independent 
and therefore potentially physical information in those at all. We shall give 
an affirmative answer by showing that (the locations of) certain singularities 
of the full propagators are indeed gauge independent Q . 
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Non-Abelian Gauge-Boson Propagator We begin by analysing the Lor- 
entz structure of the gluon propagator in the case of a general gauge that 
preserves the rotational symmetry. Moreover, we shall simplify things by 
dropping any color indices, which presupposes the absence of color symmetry 
breaking .[] 

In momentum space one can define a transverse projection of the 4- 
velocity of the heat bath, — (g^" — ^jr-)it CT , and use it to write the 
general structure of the gauge-boson propagator as 

G^(fc) = A A A^ + A B W V + A c &" + AuD^ (65) 

with 



n 

I h n k v (66) 

e^(fc) = — {n^W + Wn"} , D^(fc) = — -. 

|fc| Ar 

Here A^ 11 is the spatially transverse projector introduced already in (fl8|), 
and 23^ is a second, independent tensor that is likewise transverse with 
respect to 4-momentum, but longitudinal with respect to 3-momentum. C*" 
and complete the basis of symmetric tensors, with C M!y chosen such that 
k^Q^kv — 0, and I)'"' longitudinal with respect to 4-momentum. 

A, 23, and D are idempotent, whereas C 2 = — (23 + 23). Under a Lorentz 
trace, products of one such tensor with a different one vanish; without trace, 
A is orthogonal to all the others, but among the rest one only has 23 _L D. 

Similarly, we shall decompose the gluon self-energy II ^ = G~ lAW — Gq 1 ^ 
according to 

= -n A A^ - n B w - n c & lv - n D D^ . (67) 

At momentum scales co,k <C T, fi, the leading-order term in the one-loop 
polarisation tensor 77 Ml/ is given by the HTL (HDL) ~ max(T 2 ,/i 2 ), which 
has only 4-d-transverse contributions 

77 HT L= l (77 HTL / _ /7 HTL ) ( g ga) 
= -gijHTL (68b) 

k 

i7g TL = (68c) 
i7g TL = (68d) 



7 



The extension of the following results to color superconducting situations has 
not yet been worked ou t, bu t would be of great interest in view of the gauge 



dependence issues there [111] 
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where, at high T, 

2 T 2 

i7 HTL / = (69) 

^ 00 ~ 3 ^ 2|fe| ln p-|fc|J- (70j 

As a function of frequency and 3-momentum, the result is identical in QED 
and QCD, if for the latter we define e 2 := g 2 {N + N f /2) (for SU(/V) with N f 



flavors) [|75,124|. If there is also a nonvanishing chemical potential fi, a similar 
result holds where T 2 — > T 2 + -^i^ 2 in terms oc Nf (all of them in QED) — 
with obvious generalization to the case of different chemical potentials \Xi for 
different flavors i. 

The HTL-dressed propagator G HTL = (G 1 + /I 1 ™)" 1 has poles off the 
usual light-cone, which come in two branches determined by 

Zi™- 1 = k 2 - 7TH tl = (71a) 

ZiHTL-l = k 2 _ jjHTL = q (?lb) 

Since, as we have seen above, the HTL contribution is completely gauge 
independent and the gauge fixing parameters contained in G^ 1 do not appear 



in (71a, 71b), the A- as well as the B-part of the HTL propagator is completely 
gauge independent. 

The physical interpretation of the A- and S-branch of propagator poles 
(displayed in Fig. ||) is that the former represents quasi-particles which are 
in-medium versions of the physical polarisation of the gauge bosons, while the 
appearance of the i?-branch is a purely collective phenomenon corresponding 
to charge density oscillations (plasmons) above the plasma frequency and to 
charge screening below. 

Beyond the HTL approximation and in non-Abelian theories, however, 
one has gauge parameter dependences within 77, and also 77 M1/ fc M ^ so that 

Considering a general, rotationally invariant gauge f' J '(k)A ll (k) as in (p6[), 
this can be parameterized as 

/"(7c) = (3{k)k^ + 7(fc)n M . (72) 

Covariant gauges then correspond to = 1,7 = 0, Coulomb gauges to (3 = 
h 2 , 7 = — k°, and temporal axial gauge to (3 = k° /k 2 , 7 = 1. 

The structure functions of the gauge propagator become more compli- 
cated, to wit, 

A A = [k 2 -II A ]- 1 (73a) 

Tc -an 2 +^ 
f3 2 k 2 - aII D 



a b = [k 2 - n B - m*\no-«n* + ?f?n D] _ x (m) 



a c = J^\-^o (73c) 
p 2 k 2 - a n D y ' 
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Fig. 6. The location of the zeros of Z\^ TL_1 (transverse gluons) and of Z\^ TL_1 
(longitudinal plasmons) in quadratic scales such as to show propagating modes and 
screening phenomena on one plot. Above a common plasma frequency tv p i. there 
are propagating quasi-particle modes, which for large momenta in branch A tend to 
a mass hyperboloid with asymptotic mass = fwpj , and in branch B approach 
the light-cone exponentially with exponentially vanishing residue. For u) < w p i., 
|k| is the inverse screening length, which in the static limit vanishes for mode A 
(absence of magnetostatic screening), but reaches the Debye mass, rh 2 D = 3u)ii , for 
mode B (electrostatic screening) 



j*h* + a(k 2 -n B ) 
ZId = lib 

I3 2 k 2 - an D 



(73d) 



and there are gauge parameters everywhere, both explicitly and also within 
the structure functions of 77. 

These gauge dependences are controlled by the gauge dependence identi- 
ties discussed above, and, in compact notation, they have the form 



8A l 



.7=0 



= -(A im 5X j m + SX t m A mj ) 



j=o 



(74) 



for the full propagator. Specialized to the thermal gauge-boson propagator 
in /''-gauge, one finds |p3,p3| 



Ma = A t 



5A B l = A-/ 



fi^ 7/3-a7T c /|A; 



= A-/SY 



p 2 k 2 - an. 



D 



2n v 8X" u 



A^SZ 



(75a) 
(75b) 



but no such relations for Ac and A rj . 

If 5Y and 5Z are regular on the two "mass-shells" defined by A^ 1 = 
and Ag 1 — 0, the relations (75a, 75t) imply that the locations of these 
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particular singularities of the gluon propagator are gauge fixing independent, 
for if A ^ = = Ag 1 then also A^ 1 + SA^ 1 = = A^ + SA B \ 

So everything depends on whether the possible singularities of 6X" could 
coincide with the expectedly physical dispersion laws A^ 1 = and A^ 1 = 0. 
Diagrammatically, 5X V is obtained from the primary diagrams of Fig. ||| by 
inserting one additional vertex in all possible ways (and omitting all result- 
ing tadpole-like diagrams in the case of no spontaneous symmetry breaking). 
Since the primary diagrams are 1-particle reducible with respect to the ba- 
sic ghost line attached to SF", 8X y ^ will have singularities like the (full) 
ghost propagator. These singularities are however generically different from 
those that define the spatially transverse and longitudinal gauge-boson quasi- 
particles. Indeed, in leading-order thermal perturbation theory, the temper- 
ature power-counting rules referred to in Sect. ^| imply that the ghost prop- 
agator does not receive contributions ~ e 2 T 2 and thus will have completely 
different dispersion laws. 

The other parts of the diagrams making up 5X" are 1-p-i and may develop 
singularities for other reasons, namely when one line of such an 1-p-i diagram 
is of the same type as the external one, and the remaining ones are massless. 
This may potentially give rise to infrared or mass-shell singularities. However, 
these singularities will be absent as soon as an overall infrared cut-off is 
introduced, for example by considering first a finite volume. In every finite 
volume, this obstruction to the gauge-independence proof is then avoided, 
and A~ A X = and A^ 1 = define gaug e- ind ependent dispersion laws if the 
infinite volume limit is taken last of all [ 114 1 . 

This reasoning leads to the conclusion that the positions of all the sin- 
gularities of A a are gauge- fixing independent, though not necessarily their 
type or e.g. their residues if they are simple poles. In the case of Ap , there is 
a slight complication by the contents of the square bracket in ( 75b ). There is 
a kinematical pole 1/fc 2 hidden in the h's, and there is a contribution from 
the obviously gauge-dependent Ad (cf. (73d)). These gauge artefacts have 
to be excluded, but they are gauge dependent already at tree level and thus 
easy to identify. For example, 4^ as defined above has a factor of k 2 which 
cancels in the Coulomb gauge propagator but not in that of covariant gauges, 
so this massless mode is a gauge mode and thus unphysical. 

The gauge-(in)dependence identities (75a, 75bj ) also explain the gauge de- 
pendences found in the one- loop calculation of the plasmon damping constant 
mentioned above. Truncating e.g. (75a) at one-loop order gives 

5A A 1 ^(k) = A A 1 ^5Y^ 1 \ (76) 

with superscripts referring to bare loop order and using that SY^ = 0. 
However, the HTL plasma dispersion law is derived from A A 1(0) +4-1(1) =0 , 
and the "correction" A^ 1 ^ ~ 7T HTL ~ g 2 T 2 ~ uj 2 y is not small but sets the 
scale for everything. The temperature-power-counting rules of Sect. ^ give 
<jy(i) ~ g 2 T/ojpi., so the right-hand side of ( |76| ) does not vanish at the order 
of the damping contribution 7 x u> p i. ~ gu} 2 ^ . 
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On the other hand, if one does have a good expansion parameter (which 
bare loop order obviously is not), then the identities 
by-order gauge independence. 

As will be discussed further below, HTL perturbation theory [ |38|]39| ] claims 
to be a systematic framework, although not up to arbitrarily high orders, 
and the expansion parameter is essentially \fg^ . In |36[ ] , the long- wavelength 
plasmon damping constant has been calculated by Braaten and Pisarski with 
the result 7(|fe| = 0)/[ g 2 JJ v ] = +6.635 . . . and formal checks as to its gauge 
independence were positive. 

More explicit calculations by Baier et al., however, revealed that, in co- 
variant gauges and on plasmon-mass-shell, HTL-resummed perturbation the- 
ory still leads to explicit gauge dependent contributions to the damping of 
fermionic |l7j as well as gluonic quasi-particles. But, as was pointed out 
in |114| , these apparent gauge dependences are avoided if the quasi-particle 
mass-shell is approached with a general infrared cut-off such as finite vol- 
ume, and this cut-off lifted only in the end. This procedure defines gauge- 
independent dispersion laws and the gauge dependent parts are found to 
pertain to the residue, which at finite temperature happens to be linearly 
infrared singular in covariant gauges, rather than only logarithmically as at 
zero temperature, due to Bose enhancement. 



75a . 75b|) imply order- 



Extension to Fermions The fermion propagator at non-zero temperature 
or density has one more structure function than usually. In the ultrarela- 
tivistic limit where masses can be neglected, the fermion self-energy can be 
parametrized according to 

E(k , k) = a(k , k) 7 ° + b(k Q , k)k ■ 7 (77) 

with k = k/\k\ (again neglecting the possibility of color superconductivity). 

Defining S±{k^, k) = b(ko, k) ± a(fco, k), a natural decomposition of the 
fermion self-energy and propagator S^ 1 = — ^ + S is given by 

7o£(fco,k) = S+{k ,k)A+{k) - S-(ko,k)A-(k), (78) 

7o S- 1 (fc ,k) = A^(k 0) k)A+(k) + AzHko,k)A-(k) (79) 

with A± = -[k T (k + S±)] and spin matrices 

1 ± ■ k 

A ± (k) = ^— , A+ + A- = l, (80) 

A\ = A±, A+A- = A_A+ = 0, Tryl ± = 2, (81) 

projecting onto spinors whose chirality is equal (A + ), or opposite (^1-), to 
their helicity. 

In the HTL approximation where |fco|, |k| <ti max(T, /i), the fermion self- 
energy has been first calculated by Klimov ^(J as 
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where M 2 is the plasma frequency for fermions, i.e., the frequency of long- 
wavelength (k — > 0) fermionic excitations 



M 2 = lS± 



(83) 



(C f = (N 2 - l)/2N in SU(iV) gauge theory, and .9 2 C*/ -► e 2 in QED.) 

This leads to two separate branches of dispersion laws of fermionic quasi- 
particles A±(u), k) = carrying particle and hole quantum numbers re- 
spectively [125 107 . As shown in Fig. fj], the (— )-branch, which is occasionally 
nicknamed "plasmino" , exhibits a curious dip reminiscent of the dispersion 
law of rotons in liquid helium. 




Fig. 7. The location of the zeros of Zijjj 1 in the HTL approximation in quadratic 
scales. The additional collective modes of branch (— ) ( "plasminos" ) has a minimum 
of us at u/M « 0.93 and \k\/M S3 0.41 and approaches the light-cone for large 
momenta, but with exponentially vanishing residue. The regular branch approaches 
a mass hyperboloid (here a straight line parallel to the diagonal) with asymptotic 
mass 



As we have seen in Sect. || gauge dependences start at order T in the 
high-temperature expansion. While the HTL result (82) is completely gauge 
independent, gauge parameter dependences enter at subleading order. The 
gauge dependence identity for the fermion propagator is, however, somewhat 
simpler than that for the gluon propagator. The singularities in the fermion 
propagator can be summarized by 

det(S-i) = (84) 
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where a, g are spinor indices. The gauge dependence identity thus takes the 
form 

6 det (k) = - det (fc) [<DT r (fc) + (85) 

where the two expressions within the square brackets are Dirac traces of the 
diagrams appearing within the braces in Fig. |^. 

The conclusion of gauge independence of the solutions of ( |S4| ) can now 
be reached by essentially the same arguments as those for the gauge boson 
propagator (with the exception that now there are no additional, gauge- 
dependent kinematical poles like those arising in the projection onto mode B 



in ( 75b )). The only obstruction to gauge independence comes from singular- 
ities of 6X T T (k) and SX T f (k). If mass-shell singularities from massless gauge 
modes are avoided by infrared regularization or finite volume, only the singu- 
larities of the ghost propagator need to be considered. Again, the latter are 
generically different from those leading to the now fermionic quasi-particles 
since ghosts do not have HTL self-energies ~ T 2 . 

We have thus seen that all the singularities of the fermion propagator as 
well as those of the A- and B-branch of the gluon propagator (with some 
exceptions for the latter) are gauge-fixing independent. On the other hand, 
residues (if those singularities are simple poles at all) are not protected and 
may be gauge dependent; even the nature of the singularity may be different 
from gauge to gauge, as is well known to be the case already for the electron 
propagator in zero-temperature QED pl|. 



6 Quasiparticles in HTL Perturbation Theory 

We have already seen that loop order in bare perturbation theory is not a good 
expansion parameter for calculating corrections to quasi-particle properties 
at soft scales ~ gT. Technically what happens is that the HTL contributions 
to one-loop vertex functions are of the same order of magnitude as their 
tree-level counterparts for external momenta k ~ gT: 



HTL N^iZ-N „N-2i4-N 



Zjy ~ g T k ~g 



d N C 



dA N 



fc~gT 



Therefore all HTL contributions need to be resummed in Feynman diagrams 
that are sensitive to the soft regime k ~ gT. 

Since HTL's are the leading contributions from hard momenta k ~ T, this 
can be understood as the transition from the bare Lagrangian to an effective, 
Wilson-renormalized one for k ~ gT, C — > C + £ HTL . £ HTL is the effective 
Lagrangian containing all the HTL diagrams and arises from integrating out 
all hard modes. 

Soon after the identification of all the HTL's of QCD in |9|,|§, it has 
been found that, formally, £ HTL has a comparatively simple and manifestly 
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gauge-invariant integral representation [120,41 

£ HTL = ^2 



4tt V1 v-D{A) 



4> 



da 



-in 







(v ■ D adj .(A))- 



-F, 



(87) 



= ±77™" the 



with M 2 the fermionic plasma frequency given in ( |83| ) and < 
more familiar one of the gauge bosons (cf. (p9|)). In this integral representation 
v = (l,v) is a light-like 4-vector, i.e. with tr = 1, and its spatial components 
are averaged over by da v . v is the remnant of the hard plasma constituents' 
momenta p M ~ Tv 1 - 1 , namely their light-like 4- velocity, and the overall scale T 
has combined with the coupling constant to form the scale of thermal masses, 
M,lo p i. ~ gT. 

The HTL effective Lagrangian (|87|) is manifestly gauge invariant and 
moreover gauge independent (M and uo p i. do not depend on the gauge fix- 
ing parameters used to integrate out the hard modes). It is non-local and 
Hermitian only in a Euclidean form, i.e. prior to analytic continuation to 
real time/frequencies. It has cuts which physically correspond to the phe- 
nomenon of Landau damping. The equations of motions associated with ( pT[ ) 
can also be obtained from kinetic theory, which is extremely useful to gain 
further physical insight pj, 22], 79|,p6| . There is also a noteworthy connection 
to Chcrn-Simons theory |52|,|53||l03|| 

Using (|8^) as an effective theory for soft scales ~ gT means that the bare 
propagators are to be replace by those of HTL quasi-particles, and these have 
infinitely many nonlocal vertices. E.g., the three-gluon vertex becomes 



r°Zf +ilTL (k, q, r) = igf abc {g^(k - q) e + cycl. 



3cj 



pi. 



da v 

An 



■Vfj,V v V g 



''() 



k ■ v r ■ v k ■ v q ■ v 



In QCD, there are HTL vertices for any number of external gluons and up to 
two quark lines, whereas in QED, where v ■ D a dj.(A)) — > v ■ d in (|S7|), there is 
"only" an HTL photon self-energy 77 M „, an HTL fermion self-energy E, and 
vertices involving two fermions and an arbitrary number of photons. 

While the effective Lagrangian (|87j ) is gauge invariant and gauge inde- 
pendent in its entirety, NLO corrections won't be so. However, as we have 
seen in the previous section, the positions of singularities of the effective 
(quasi-particle) p ropa gators are protected against gauge dependences by the 
identities (^5aj), and (ph. 
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6.1 Long- Wavelength Plasmon Damping 

The first such correction to be calculated by means of the HTL-resummed 
perturbation theory was the damping rate of long- wavelength plasmonsn from 
the shift of the pole of the gluon propagator at k = from u> = ujpf L — > 
uj = LUp\, — i7(fc = 0) with the result |36| 

7 (fc = 0) = +6.635 . . . ^-JP- = 0.264%/]% ^ TL , (89) 

247T 

implying the existence of weakly damped plasmons for jCl. 

In QCD, where one is interested in the range g ~ 1, one finds that the ex- 
istence of plasmons as quasi-particles requires that g is significantly less than 
2.2, so real QCD is on the borderline of having identifiable long-wavelength 
quasi-particles. 

The corresponding quantity for fermionic quasi-particles has been calcula- 
ted in J8l],[40) with a comparable result: weakly damped long- wavelength 
fermionic quasi-particles in 2- or 3-flavor QCD require that g is significantly 
less than 2.7. 



6.2 NLO Correction to Gluonic Plasma Frequency 



In [117 1, Schulz has calculated also the real part of the NLO contribution 
to the gluon polarization tensor in the limit of k — *■ which determines the 
NLO correction to the gluonic plasma frequency. 

The original power-counting arguments of [|39| suggested that besides one- 
loop diagrams with HTL-resummed propagators and vertices, there could be 
also contributions from two- loop diagrams to relative order g. The explicit 



(and lengthy) calculation of [117| showed that those contribute only at order 
g 2 ln(l/g) rather than g, and the NLO plasma frequency in a pure-glue plasma 
was obtained as 



, — , , HTL 



1 - 0.09V Ng 



(90) 



In this particular result, HTL-resummed perturbation theory turns out 
to give a moderate correction to the leading-order HTL value even for g ~ 1; 
see however below. 

While the calculations leading to ( p9| ) and (^0|) contain some interesting 
physics, in the following we shall go into more detail only for a couple of 
more tractable cases, which nonetheless will turn out to involve a number of 
salient points. 

8 For kinematical reasons, there should be no difference between spatially trans- 
verse and longitudinal gluonic quasi-particles (cf. Fig. |J), since with k — > one 
can no longer tell the one from the other. However, the limit k — *■ involves 
infrared problems (see further below), and there are even explicit calculations 
that claim to find obstructions to this equality, which are however refuted 
by the recent work of S . 
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6.3 NLO Correction to the Non-Abelian Debye Mass 

The poles of the (gluon) propagator do not only give the dispersion law of 
quasi-particles, but also the screening of fields with frequencies below the 
plasma frequency and in particular of static fields. Below the plasma fre- 
quency, there are poles for k 2 < 0, as displayed in Fig. ^, corresponding in 
configuration space to exponential fall-off with (frequency-dependent) screen- 
ing mass •v/|fe 2 ], i.e. screening length l/-v/|fe 2 ]. 

In the static case, branch A of the gluon propagator describes the screen- 
ing of (chromo-)magnetostatic fields. While there is a finite screening length 
as long as a; > 0, the A-branch of the HTL propagator becomes unscreened in 
the static limit. Whereas in QED, a "magnetic mass" is forbidden by gauge 
invariance (5^,|2^], some sort of entirely non-perturbative magnetic mass is 
expected in non-Abelian gauge theories in view of severe infrared problems 



caused by the self-interactions of magnetostatic gluons |110, 96|j64|| . 

Branch B 1 on the other hand, contains the information about screening of 
(chromo-)electric fields as generated by static charges (Debye screening). The 
Debye mass given by the leading-order HTL propagator is ihjj — \/3wpi.- The 
determination of its NLO correction has a history that is at least as long as 
the plasmon (damping) puzzle, for it starts already with (ultra-relativistic) 
QED. 

Customarily, the Debye mass (squared) has been defined as the infrared 
limit 7Too(w = 0, k — > 0), which indeed is correct at the HTL level, cf. ( |68b[ ) 
and @. 

In QED, this definition has the advantage of being directly related to 
a derivative of the thermodynamic pressure, so that the higher-order terms 
known from the latter determine those of ilgo (uj = 0, k — ► 0) through 
58|j7e 



d 2 P 



m=o d/i 2 



e 2 T 2 ( 3e 2 y/% 



3 



77 00 (0,fc^0) =e 2 — =—- \ — + + (91 



M =o 3 \ 8tt 2 4tt 3 



This result is gauge independent because in QED all of is. 

In the case of QCD, there is no such relation. In fact, one expects 8m 2 D /rh 2 D 
~ g rather than g 3 because of gluonic self-interactions and Bose enhancement. 
The calculation of this quantity should be much easier than the dynamic 
ones considered above, because in the static limit the HTL effective action 
collapses to just the local, bilinear HTL Debye mass term, 

£ HTL static _l^ tlvl 2 (g2) 

This is also gauge invariant, because Aq behaves like an adjoint scalar un- 
der time-independent gauge transformations. Resummcd perturbation theory 
thus boils down to a resummation of the HTL Debye mass in the electrostatic 
propagator, which is what had been done already since long [|63|j7q|. 
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Using this simple ("ring") resummation in QCD, one finds however the 



gauge dependent result [121 




n O0 (0,0)/mj D = l + a—J———g (93) 



where a is the gauge parameter of general covariant gauge (which coincides 
with general Coulomb gauge in the static limit). 

This result was interpreted as meaning that the non-Abelian Debye mass 



could not be obtained in resummed perturbation theory [ 1 1 1 or that one 
should use a physical gauge instead [ [73|]7^] . In particular, temporal axial 
gauge was put forward, because in this gauge there is, like in QED, a linear 
relationship between electric field strength correlators and the gauge propaga- 
tor. However, because static ring resummation clashes with temporal gauge, 
inconclusive and contradicting results were obtained by different authors 
[|72U61|J73|, and in fact one cannot do without vertex resummations if one 



wants to be consistent there 15,106]. But be that as it may be, switching 



to the chromoelectric field strength correlator is not good eno ugh, for it is 



gauge variant and its infrared limit is equally gauge dependent |116 |. 

On the other hand, in view of the gauge dependence identities discussed 
in the previous section, the gauge dependence of (|93| ) is no longer surprising. 
Gauge independence can only be expected "on-shell" , which here means u> = 
but k 2 -rh%. 

Indeed, the exponential fall-off of the electrostatic propagator is deter- 
mined by the position of the singularity of Z\_b(0, k), and not simply by its 
infrared limit. This implies in particular that one should use a different def- 
inition of the Debye mass already in QED, despite the gauge independence 
of (fl]), namely gl| 

mf> = 7T o(0,fc) (94) 

k — > — 

For QED (with massless electrons) , the Debye mass is thus not given by 
([»]) but rather as 

m 2 D = i7 00 (0, k -> 0) + [7T 00 (0, fc)| fe2= _ mi - #oo(0, k - 0)] 

e 2 T 2 / 3e 2 V3e 3 e 2 .n 4. \ . R . 

where jl is the renormalization scale of the momentum subtraction scheme,^ 



i.e. n^(k 2 = -A 2 )| T =o = 0. Since de/dln/2 = e 3 /(127r 2 ) + 0(e 5 ), (1951) is a 



renormalization-group invariant result for the Debye mass in hot QED, which 
( |9ll ) obviously failed to be. 

9 The slightly different numbers in the terms oc e 4 T 2 quoted in p^ , p3[ ] pertain to 
the minimal subtraction (MS) scheme. 
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In QCD, where gauge independence is not automatic, the dependence on 
the gauge fixing parameter a is another indication that (|9^) is the wrong 
definition. For (Q) we need the full momentum dependence of the correction 
<577 o(fco = 0, k) to i7^, TL . Since only the electrostatic mode needs to be 



dressed, this is not difficult to obtain [115] 

sn 00 



i f <P~ 2e p 

)(fco = 0, k) = gm D N ^ ^ivflvj J (27r )3-2 £ 



9 2 T 

1 1 4rh 2 D - (k 2 + m 2 D )[3 + 2pk/p 2 



p 2 + m 2 D p 2 P 2 [(p + k) 2 + m 2 D ] 

+^ 2 +^) p4[( ;^ ) y + fc ^ ] }- (^) 

In accordance with the gauge dependence identities, the last term shows 
that gauge independence holds algebraically for k 2 = —rh 2 D . On the other 
hand, on this "screening mass shell" , where the denominator term \(p + 
k) 2 + m 2 D \ — > [p 2 + 2pk], we encounter IR-singularities. In the a-dependent 
term, they are such that they produce a divergent factor 1/ [k 2 + m 2 D ] so that 
the gauge dependences no longer disappear even on-shell. This is, however, 
the very same problem that had to be solved in the above case of the plas- 
mon damping in covariant gauges. Introducing a temporary infrared cut-off 
(e.g., finite volume), does not modify the factor [fc 2 + m 2 D \ in the numerator 
but defuses the dangerous denominator. Gauge independence thus holds for 
all values of this cut-off, which can be sent to zero in the end. The gauge 
dependences are thereby identified as belonging to the (infrared divergent) 
residue. 

The third term in the curly brackets, however, remains logarithmically 
singular on-shell as the infrared cut-off is to be removed. In contrast to the in- 
dependent term, closer inspection reveals that these singularities are coming 
from the massless magnetostatic modes and not from unphysical massless 
gauge modes. 

At HTL level, there is no (chromo-)magnetostatic screening, but, as we 
have mentioned, one expects some sort of such screening to be generated 
non-perturbatively in the static sector of hot QCD at the scale g 2 T ~ gmo 



While this singularity prevents evaluating (|96|) in full, the fact that this 
singularity is only logarithmic allows one to extract the leading term of (Bq) 



under the assumption of an effective cut-off at p ~ g T as HE] 



Sm 2 D 



D 




0(g). (97) 



The 0((7)-contribution, however, is sensitive to the physics of the magne- 
tostatic sector at scale g 2 T, and is completely non-perturbative in that all 
loop order > 2 are expected to contribute with equal importance. 
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Because of the undetermined 0(<?)-term in (p7\), one-loop resummed per- 
turbation theory only says that for sufficiently small g, where 0(g\n(l/ g)) 3> 
O(g), there is a positive correction to the Debye mass of lowest-order pertur- 
bation theory following from the pole definition (|94|), and that it is gauge 
independent. 

On the lattice, the static gluon propagator of pure SU(2) gauge theory at 
high temperature has been studied in various gauges fl68p(| with the result 
that the electrostatic propagator is exponentially screened with a screening 
mass that indeed appears to be gauge independent and which is about 60% 
larger than the leading-order Debye mass for temperatures T/T c up to about 
10 4 . 



In [116], an estimate of the 0(g) contribution to (97) has been made using 
the crude approximation of a simple massive propagator for the magnetostatic 
one, which leads to 



4 



D 




2m D 1 
In - 

m m 2 



(98) 



On the lattice one finds strong gauge dependences of the magnetostatic 
screening function, but the data are consistent with an over-all exponen- 
tial behaviour corresponding to m m w 0.5g 2 T in all gauges |6^ , ^5| . Using 
this number in a self-consistent evaluation of (pq) gives an estimate for mjj 
which is about 20% larger than the leading-order value for T/T c = 10 . . . 10 4 . 

This shows that there are strong non-perturbative contributions to the 
Debye screening mass mu even at very high temperatures. Assuming that 
these are predominantly of order g 2 T, one-loop resummed perturbation the- 
ory (which is as far as one can get) is able to account for about 1/3 of this 
inherently non-perturbative physics already, if one introduces a simple, purely 
phenomenological magnetic screening mass. 



Other Non-Perturbative Definitions of the Debye Mass A differ- 
ent approach to studying Debye screening non-perturbatively without the 
complication of gauge fixing is to consider spatial correlation functions of 
appropriate gauge-invariant operators such as those of the Polyakov loop 



L(x) = -^IrPexp j-uj d T A (T,x) 



(99) 



The correlation of two such operators is related to the free energy of a quark- 
antiquark pair [ 100 1 . In lowest order perturbation theory this is given by the 
square of a Yukawa potential with screening mass rhu [ 101 1 ; at one-loop order 
one can in fact identify contributions of the form ( p8[ ) if one assumes magnetic 
screening p3,116|, but there is the problem that through higher loop orders 
the large-distance behaviour becomes dominated by the magnetostatic modes 
and their lightest bound states p3 . 
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In |I(J, Arnold and Yaffe have proposed to use Euclidean time reflection 
symmetry to distinguish electric and magnetic contributions to screening, 
and have given a prescription to compute the sublogarithmic contribution 
of order g 2 T to mjj nonperturbatively. This has been carried out in 3-d 
lattice simulations for SU(2) as well as for SU(3) @. The Debye mass 

thus defined shows even larger deviations from the lowest-order perturbative 
results than that from gauge-fixed lattice propagators. E.g., in SU(2) at T — 
10 4 T C this deviation turns out to be over 100%, while in SU(3) the dominance 
of g 2 T contributions is even more pronounced. 

Clearly, (resummed) perturbation theory is of no use here for any tem- 
perature of practical interest. However, the magnitude of the contributions 
from the completely nonperturbative magnetostatic sector depends strongly 
on the quantity considered. It is significantly smaller in the definition of the 
Debye mass through the exponential decay of gauge-fixed gluon propagators, 
which, as we have seen, leads to smaller screening masses on the lattice (and 
gauge-independent ones, too, apparently). In quantities where the barrier in 
perturbation theory arising from the magnetostatic sector occurs at higher 
orders, HTL-resummed perturbation theory should be in much better shape, 
and we shall find some support for this further below. 



6.4 Dynamical Damping and Screening 



A logarithmic sensitivity to the nonperturbative physics of the magnetostatic 
sector has in fact been encountered early on also in the calculation of damping 
of a heavy fermion [ 108 , and more generally of hard particles Q,^5 42,112]. It 
also turns out to occur for soft quasi-particles as soon as they are propagating 
and not just stationary plasma oscillations. 



Because this logarithmic sensitivity arises only if one internal line of (re- 
summed) one-loop diagrams is static, the coefficient of the resulting g ln(l/p)- 
term is almost as easy to obtain as in the case of the Debye mass, even though 
the external line is non-static and soft, requiring HTL-resummed vertices (see 
Fig. 1). 




1/2 "^^^^^ + 1/2 

Fig. 8. One-loop diagrams in HTL-resummed perturbation theory. HTL-resummed 
quantities are marked with a blob 



The infrared singularity arises (again) from the dressed one-loop diagram 
with two propagators, one of which is magnetostatic and thus massless in 
the HTL approximation, and the other of the same type as the external one, 
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so only the first diagram in Fig. g is relevant. The dressed 3- vertices in it 
are needed only in the limit of one leg being magnetostatic and having zero 
momentum. Because of the gauge invariance of HTL's, these are determined 
by the HTL self-energies through a differential Ward identity, e.g. 



d 

r^vgik; — fc; 0) = — -^^n ^(k) 



for the 3-gluon vertex (color indices omitted). 
Comparatively simple algebra gives p6| 

Sn T (k) ~ -g 2 N4k 2 [l + d k 2lI I (k)} 2 S I (k), 

where 

d 3 p 1 -1 



(2tt) 3 p 2 {k-p) 2 - n I (k-p) 



A.B 



k 2 =n I (k), P °=o 



(100) 



(101) 



(102) 



and the logarithmic (mass-shell) singularity arises because (k — p) 2 — IIj(k - 



P) 



-> -p 2 + 2pk - 77/(fc -p) + 77/(fc) ~ \p\ as k 2 -> ilj(fe). 
The IR-singular part of S/(fc) is given by 

h ( k) = t r ** 1 



(2tt) 3 p 2 p 2 - 2pk + IIi(k-p)- III (k) - ie 
d 3 p 1 1 
(27r) 3 p 2 p 2 — 2pk — ie 
50 dp 1 p + 2|fc|-is 
~p2\k\ %- 2\k\ - ie 



Tll + ^TT^fc)]- 1 
Tll + d^IIiik)]- 1 J 



(103) 



where in the last line we have inserted an IR cutoff A <C gT for the p-integral 
in order to isolate the singular behaviour. 

One finds that (103) has a singular imaginary part for propagating modes, 



Sz(*) c + d^n^k)}- 1 In M + O(A ) 

B7T fc A 



(104) 



for A; 2 > (fc real), and a singular real part in screening situations where 
|fc| 2 — > — k 2 , k£1 (i.e., fc imaginary): 



S/(fc) ~ + dk.n^k)]- 1 In £ + O(A ). 



(105) 



So from one and the same expression we can see that logarithmic IR 
singularities arise whenever |fc| ^ 0, leading to IR singular contributions to 
damping or (dynamical) screening, depending on wheth er uj > or < and 
thus fc 2 > or < 0. The case fc = is IR-safe, because (101) is proportional 
to fc 2 , while 



S/(fc) 



T /Tlfcl 

[i + d^n^k)]- 1 + o' 



4tt 2 A 



A 2 



for fc -► 0. 



(106) 
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This shows that there is a common origin for the infrared sensitivity of 
screening and damping of HTL quasi-particles. The perturbatively calculable 
coefficients of the resulting g ln(l/p)-terms are in fact beautifully simple: For 



the damping of moving quasi-particles one obtains [109 56 
9 2 NT 



7/(1*1) 



Air u(\k\)[l-d t 



dyn^k)] 1 = g NT 
iTi(fc)] g~ 4tt 



u/(|fe|)ln- 



9 



(107) 
= 0). 



where uj(|fe|) is the group velocity of mode / (which vanishes at k 
The IR-sensitive NLO correction to screening takes its simplest form when 
formulated as flsll 



9 2 NT 
2tt 



1 



k 7 (w) In - +0(1) 



9 



(108) 



where Ki(u>) is the inverse screening length of mode / at frequency u> < tv p i, 
(which in the static limit approaches the Debye mass and perturbatively 
vanishing magnetic mass, resp., while approaching zero for both modes as 

UJ — > LOpl.). 

A completely analogous calculation for the fermionic modes (for which 
there are no screening masses) gives 



7±(|fc|) = 



g 2 C F T t 

4"7T 



\k\ 



In - + 0(1) 
9 



(109) 



for \k\ > 0. The group velocity v± equals ±| in the limit (|fe|) — * 0, and 
increases monotonically towards +1 for larger momenta (with a zero for the 
(-)-branch at \k\/M ps 0.41). For strictly |fe| = 0, the IR sensitivity in 
fact disappears because (109) is no longer valid for |fc| -C A, but one has 
7±(|fc|)|sing. oc g 2 T\k\/X instead. Thus 7±(0) is calculable at order g 2 T in 
HTL-resummed perturbation theory, and has been calculated in |81 40||. 



The fermionic result (109) applies in fact equally to QED, for which one 
just needs to replace g 2 Cj? — > e 2 . This is particularly disturbing as QED does 
not allow a non-zero magnetic mass as IR cutoff, and it has been conjectured 
that t he d amping 7 ~ g 2 T or e 2 T itself might act as an e ffective IR cutoff 
[ p4|j95 ,109,4|, which however led to further difficulties [ 105 1 . The solution for 
QED was finally found by Blaizot and Iancu [f^p4p5|| who showed that 
there the fermionic modes undergo over-exponential damping in the form 
e -7* _> e -^Ttin(u ph t) ( for v 1^ so g nite time j g the actua i cut-off. The 
fermion propagator has in fact no simple quasi-particle pole, but nevertheless 
a sharply peaked spectral density. 

In non-Abelian theories, on the other hand, one does expect static (chro- 
mo-)magnetic field to have finite range, and lattice results do confirm this, 
so the above estimates may be appropriate after all, at least for sufficiently 
weak coupling. 



Thermal Gauge Field Theories 



37 



6.5 NLO Corrections to Real Parts of Dispersion Laws 

The above analysis has identified the imaginary parts of the dispersion laws 
to be sensitive to non-perturbative IR physics except at k = and where 
the group velocity vanishes (which includes one further point at |fe| =/= 
for the fermionic plasmino branch). On the other hand, the real parts of 
the dispersion laws of gluonic and fermionic quasi-particles are IR-safe in 
NLO HTL-resummed perturbation theory. However, such calculations are 
tremendously difficult, and only some partial results exist so far in QCD 




In the following, we shall restrict our attention to the case k 2 /ui^ ^> 1 and 
consider the two branches of the gluon/photon propagator in turn. In both 
cases, interesting physics will be seen to be contained in the NLO corrections. 



Longitudinal Plasmons For momenta k 2 ^> w 2 j , the longitudinal plas- 
mon branch a pproa ches the light-cone, as can be seen in Fig. [| From k 2 — 
7Tf TL (fc) and (js8b|) one finds 

w|(|fc|) ^k 2 (l + 4fc 2 e- 6fc2 /( e2T2 >) (110) 



with e 2 = g 2 (N + Nf/2) in QCD, so the light-cone is approached exponen- 
tially. If one also calculates the residue, one finds that this goes to zero at 
the same time, and exponentially so, too. 

Instead of QCD, we shall consider the analytically tractable case of mass- 
less scalar electrodynamics as a simple toy model with at least some simi- 
larities to the vastly more complicated QCD case in that in both theories 
there are bosonic self-interactions. There are however no HTL vertices in 
scalar electrodynamics, which makes it possible to do complete momentum- 
dependent NLO calculations [ jS5| . 

Comparing HTL values of and NLO corrections to 77b , one finds that as 
k 2 — > there are collinear singularities in both: 



i7f TL (fc)/fc 2 ^^.ln^ (111) 



3 2 i ^ 2 
diverges logarithmically]^], whereas 



sn B /k 2 -+- e/i 2 c tt ±J= ( ii2) 



(with /i S c.th. oc eT the thermal mass of the scalar). Because (112) diverges 
stronger than logarithmically, one has 811 b > -ff|j TL eventually as k 2 — > 0. 
Clearly, this leads to a breakdown of perturbation theory in the immediate 
neighbourhood of the light-cone (k 2 /\k\ 2 < (e/ In \) 2 ) 1 which this time is not 

10 This is in fact the technical reason why the longitudinal branch approaches the 
light-cone exponentially when k 2 ^> . 
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caused by the massless magnetostatic modes, but rather by the massless hard 
modes contained in the HTL's. 

However, a self-consistent gap equation for the scalar thermal mass implies 
that also the hard scalar modes have a thermal mass ~ eT. Including this by 
extending the resummation of the scalar thermal mass to hard internal lines 
renders 77b regular up to and including the light-cone one obtains 



rrrcsum. 

lim B 



k 2 3fc 2 



2T 1 , C'(2) 

In h- - 7_e 



Mscth, 2 ,a C(2) 
b 



(113) 



The finiteness of (113) makes it possible that there is now a solution to the 



dispersion law with k 2 = at k 2 /(e 2 T 2 ) = |ln 2 - 09 e 4 - + 0(e). Because all 
collinear singularities have disappeared, continuity implies that there are also 
solutions for space-like momenta k 2 < 0, so the longitudinal plasmon branch 
pierces the light-cone, having group velocity v < 1 throughout, though, as 
shown in Fig. [|. While at HTL level, the strong Landau damping at k 2 < 
switches on discontinuously, it now does so smoothly through an extra factor 
exp[— e-y/|fe|/[8(|fe| — w)]], removing the longitudinal plasmons through over- 
damping for (|fc| — w)/|fe| > e 2 . 

So the collinear singularities of HTL-resummed perturbation theory on 
the light-cone were associated with a slight but nevertheless qualitative change 
of the spectrum of longitudinal plasmons: instead of being time-like through- 
out and existing for higher momenta, albeit with exponentially small and 
decreasing residue and effective mass, they become space-like at a particular 
point |fe| ~ eTTn - and expire through Landau damping soon thereafter. 

This phenomenon is in fact known to occur in non-ultrarelativistic (T < 
m e ) QED [ 122 1 , and has been considered in the case of QCD in a little-known 



paper by Silin and Ursov |119[ , who speculated that it may lead to Cherenkov 
phenomena in the quark-gluon plasma. 

In QCD, the situation is in fact much more complicated. Under the as- 
sumption that the collinear singularities are removed solely by the resum- 
mation of asymptotic gluonic and fermionic thermal masses in hard internal 
lines, the value of |fe| where longitudinal plasmons turn space- like has been 
calculated in Q. For a pure-glue plasma, it reads 

1 48 

fc 2 rit = 5 2 T 2 [ln-— - + 0(g)}. (114) 

Such an extended resummation can in fact be related to an improved and 
still gauge- invariant version of the HTL effective action |5j| , however it may 
well be that damping effects are of equal importance here (in contrast to 



scalar electrodynamics), so that (114) may not be complete. A similar un- 



solved problem occurs in the calculation of the production rate of real, non- 
thermalized photons in a quark-gluon plasma from HTL-resummed pertur- 
bation theory pLpUTl. 
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Fig. 9. The longitudinal plasmon branch of scalar electrodynamics including NLO 
corrections to the HTL result. The upper of the four lines gives the HTL result and 
the lines below correspond to NLO corrections with e = 0.3, 1, and 2, respectively. 
The latter three lines cross the light-cone such that the phase velocity starts to 
exceed 1, but with group velocity < 1 throughout. In the space-like region, the 
plasmon modes are damped by Landau damping, which is strong except in the 
immediate neighborhood of the light-cone, where it is suppressed by a factor of 
exp{-eVfe/[8(fc-w)]} 



Taken at face value, (114) would imply that propagating longitudinal 
plasmons do no longer exist for g > 1.48, and a negative O(g) contribution 
would even lower this bound. 



Energetic Quarks and Transverse Gluons and Their Role in Self- 
Consistent Thermodynamics At high momenta k 2 /ui^ ^> 1 the addi- 
tional collective modes of longitudinal plasmons and "plasminos" disappear. 
At HTL level, they do so because the residues of the corresponding poles in 
the gluon and quark propagators vanish exponentially, whereas at NLO, as we 
have just seen, they cross the light-cone and die from strong Landau damp- 
ing. The remaining transverse gluonic and normal quark modes on the other 
hand approach asymptotic mass hyperboloids. For transverse gauge bosons 
the asymptotic thermal photon/gluon mass of the HTL approximation reads 

< TL --L = ^ (H5) 

(e 2 = (N + Nf/2)g 2 for gluons); whereas in the case of fermions we have 
2\k\Sf Th -> 2M 2 (116) 



with M the HTL fermionic plasma frequency given by (J83|). 
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These results remain the correct LO ones even for u>, \k\ ~ T, because the 
light-cone values of II a and S+ are identical to their HTL /HDL values there 
and do not depend on the HTL approximation that u>, \k\ <C T ]8^ , p5[ . 

The asymptotic thermal masses play an interesting role in self-consistent 
(approximations to) thermodynamics ]2g| , |29| , |30[ : The LO (oc g 2 ) interaction 
piece of the entropy density can be expressed in terms of the light-cone values 
of the various self-energies and thus the asymptotic thermal masses. E.g., in 
the pure-glue case, the (^-contribution to the entropy density reads 



•s 



= _(^^r = -^(^i) 9 V. (nr, 

Fcrmionic contributions give similarly 

sf=-^lMlT, (118) 

possibly with nonzero chemical potential fi. With nonzero /z, one can also 
consider the quark density, which likewise is determined by the asymptotic 
mass: 

nf = -f^M^. (119) 

Up to a T- and /i- independent integration constant, entropy and quark den- 
sities determine the complete thermodynamical potential, and the above for- 
mula give nice, universal formulae for the LO interaction terms. 

Remarkably, also the NLO interaction term oc g 3 can be directly related 
to the properties of HTL/HDL quasiparticles. The so-called plasmon term of 
the thermodynamic potential oc g 3 is usually understood as arising from the 
resummation of the static Debye mass, which needs to be kept only in the 
zero modes of the electrostatic gluon propagator. The resulting coefficient of 
the order-g 3 contribution to the thermodynamic potential turns out, however, 
to have an uncomfortably large valuej^] and appears to spoil completely the 
convergence of perturbation theory for all temperatures smaller than some 
10 5 T C . 

While it is correct that all that is needed for a calculation of the ther- 
modynamic potential through order g 3 is to approximate quarks and gluons 
by their vacuum spectral densities except for the one massive electrostatic 
mode , this is clearly a cruder approximation than that of HTL-resummed 
propagators which contain a lot of physics beyond Debye screening. 

In |28|]29| , |30t it has been shown recently that in a self-consistent formu- 
lation of the thermodynamic potentials entropy and density one can find a 

The same hold s tru e for the order-p 5 contribution which has been calculated for 



QCD in UJl2,12£,3 
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real-time description of those using quasi-particles which at soft momenta 
are described by the HTL effective propagators and at hard momenta by 
their light-cone limits and NLO corrections thereof. Doing so, it turns out 
that a larger part (up to |) of the (soft) plasmon effect oc g 3 comes from 
the NLO corrections to the hard asymptotic masses, reflecting a massive^ 
reorganization of usual (Debye-screened) perturbation theory: 

=-< w2 -Vl?i^^« (120) 

(in the case of pure glue). 

6ml in HTL-resummed perturbation theory is a non-local (momentum- 
depedent) correction, which is infrared safe and thus calculable. Through the 
relation ( |120| ) one can define the average correction 

6ml = -^g 2 NTm D (121) 

which has a remarkably simple form. Similarly, for fermions one finds 

SMl = -^-g 2 C f Tm D . (122) 

Now, numerically, this correction is uncomfortably large: 



1 - g « 1 - g (123) 



(pure glue) so that perturation theory seems to become completely useless 

f° r 9 ^ 1) i- e -; a s > 0.1. 

However, a very similar problem arises already in simple scalar 4 theory. 
If one considers the large- ./V limit of the iso-vector 0(N) g 2 4> 4 theory, one 
can write down an exact gap equation of the form | 



, , d 3 k n(\/k 2 + to 2 ) 3to 2 / m 2 , 

to, = 12g A / - — — y — '- + In — - 1 (124) 

!J J (2tt)3 Vfc 2 + m 2 4tt 2 V ) 

whose solution has a perturbative expansion beginning as 

3 



2 2t-i2 



g 2 T 2 (l--g +...), (125) 



7T 



which happens to have the same O(g) coefficent as the QCD result (|12l| ), and 
which likewise gives nonsense such as tachyonic thermal masses for g>l. 



Pun intended. 
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However, if one instead writes down an approximate gap equation by 



expanding in powers of m/T and dropping terms of order (m/T) 2 



9 



9 2 T 2 - -g 2 Tm, 

7T 



(126) 



then one finds that the solution to this simple quadratic equation in m gives 
a function m(g) that is perturbatively equivalent to (125), but does not go 
mad for g > 1. On the contrary, for the standard choice of renormalization 
scale /2 = 2itT in MS, it gives a remarkably accurate approximation of the 
solution to the full gap equation (124), as is shown in Fig. filx 




Fig. 10. Various approximations to the thermal mass of a scalar boson in large- 
TV^ 4 theory: leading-order HTL (LO), next-to-leading order (NL O) a s given by 
125), and the approximately self-consistent (ASC) gap equation (126), which is 
perturbatively equivalent to NLO. The MS renormalization scale is varied by a 
factor of 2 about fi = 2nT 



Implementing analogous "approximately self-consistent" gap equations 
for the hard modes, a non-perturbative, UV finite and gauge-invariant ap- 
proximation to entropy and density of hot QCD has been proposed in [p8|j29| , |30|| 
It is perturbatively equivalent to conventional Debye-resummed perturbation 
theory but goes beyond the latter in incorporating all of the collective phe- 
nomena contained in HTL propagators as well as NLO effects in their asymp- 
totic masses. When compared to available lattice data |3^] (see Fig. [y]for the 
pure-glue case), remarkable agreement is found down to temperatures ~ 3T C . 
By contrast, conventionally resummed perturbation theory at order g 3 leads 
to S/Ssb > 1 for all but exceedingly high temperatures. 

An optimistic conclusion one could draw from this is that the transition to 
gluonic and quark quasi-particles is able to absorb a large part of the strong 
elementary interaction into the spectral properties of the former, and that, 




Fig. 11. Comparison of results from approximately self-consistent thermodynamics 
based on the H TL p ropagator and a next-to-leading approximat ion (NLA) using 
the analogue of (126) for the asymptotic gluon mass correction (121) with lattice 
data from |Q . The gray band gives roughly the lattice data with their errors. The 
analytical results are given with two boundaries corresponding to a variation of the 
MS renormalization scale p, from ttT to AnT 



at least in infrared-safe situations, these quasi-particles have comparatively 
weak residual interactions even in QCD at temperatures a few times above 
the transition temperature. 

Let me recall that even in the infrared-unsafe case of NLO corrections 
to the Debye mass, the self-consistent NLO result (|9^) using a phenomeno- 
logical magnetic mass gives the qualitatively correct result of a substantially 
increased electric mass, while underestimating the magnitude of the increase 
by a factor of 3 when compared to lattice simulations of chromoclcctrostatic 
propagators. 



7 Conclusions 

Let us summarize the findings that are of specific interest to a perturba- 
tive formulation of non-Abelian gauge theories at finite temperature and/or 
density: 

The leading-order results for self-energies and vertices in a high-temper- 
ature/density expansion, the so-called hard thermal (dense) loops, form a 
gauge-invariant and gauge-independent effective action, which is the basis for 
a systematic perturbative expansion in powers of g (rather than g 2 ), as long 
as one does not run into the perturbative barrier formed by the completely 
non-perturbative self-interacting chromomagnetostatic modes. 

Beyond leading order, gauge dependences appear in all Green functions 
of the fundamental fields. In particular the propagators which are expected 
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to carry information on gluonic or fermionic quasi-particlcs depend on gauge- 
fixing parameters. The gauge dependence identities that we have discussed 
above imply, however, that under certain conditions the location of the sin- 
gularities which define the dispersion laws of these quasi-particles are gauge- 
independent, though not e.g. residues or even the type of the singularities, 
which need not be simple poles. 

Already at NLO, screening lengths and damping constants are logarith- 
mically infrared sensitive to the nonperturbative magnetostatic sector, with 
the exception of zero 3-momentum. Infrared-safe quantities are also the real 
corrections to the dispersion laws, which in the case of longitudinal plasmons 
(and also of the plasmino branch of fermions) lead to a finite 3-momentum 
range, and, at its upper end, to space-like phase velocities. In the case of 
transverse gluonic quasiparticles and the normal branch of fermionic ones, 
the NLO corrections play an important role in self-consistent formulations of 
thermodynamics (the equation of state). 

Particularly in QCD, one faces the problem that corrections to LO results 
are rather large for almost all values of the coupling of interest. However, we 
have seen indications that this poor convergence of thermal perturbation 
theory may be overcome in approximately self-consistent reformulations .p| 
Where those can be implemented, the picture of weakly interacting quasi- 
particles even in strong interactions seems to have some support from com- 
parison with lattice data (where the latter are available), and may remain 
valid down to a few times the deconfinement phase transition temperature. 
This picture, being primarily set up in real Minkowski space, is complemen- 
tary to lattice or dimensional reduction formulations, and allows (analytical) 
calculations from first principles also where lattice gauge theory calculations 
are not (yet) feasible. Its potentials, in particular when combined with re- 
sults from other nonperturbative approaches, are, in my opinion, not yet fully 
explored. 
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